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Abstract 

We establish several fundamental properties of analysis-suitable T- 
splines which are important for design and analysis. First, we character- 
ize T-spline spaces and prove that the space of smooth bicubic polynomi- 
als, defined over the extended T-mesh of an analysis-suitable T-spline, is 
contained in the corresponding analysis-suitable T-spline space. This is 
accomplished through the theory of perturbed analysis-suitable T-spline 
spaces and a simple topological dimension formula. Second, we establish 
the theory of analysis-suitable local refinement and describe the condi- 
tions under which two analysis-suitable T-spline spaces are nested. Last, 
we demonstrate that these results can be used to establish basic approxi- 
mation results which are critical for analysis. 

Keywords: T-splines; isogeometric analysis; local refinement; analysis- 
suitable; approximation. 



1 Introduction 

T-splines were originally introduced in Computer Aided Design (CAD) as a 
superior alternative to NURBS [I] and have since emerged as an important 
technology across several disciplines including industrial, architectural, and en- 
gineering design, manufacturing, and engineering analysis. T-splines can model 
complicated designs as a single, watertight geometry and can be locally re- 
fined [U [3] . These basic properties make it possible to merge multiple NURBS 
patches into a single T-spline [H [I] and any trimmed NURBS model can be 
represented as a watertight T-spline [5] . 



The use of T-splines as a basis for isogeometric analysis has gained widespread 
attention [71 El El El HH El H21 US] • Isogeometric analysis was introduced 
in [14] and described in detail in [15]. The isogeometric paradigm is simple: 
use the smooth spline basis that defines the geometry as the basis for analy- 
sis. Traditional design-through-analysis procedures such as geometry clean-up, 
defeaturing, and mesh generation are simplified or eliminated entirely. Addi- 
tionally, the higher-order smoothness provides substantial gains to analysis in 
terms of accuracy and robustness of finite element solutions [El [TTJ [18] . 

An important development in the evolution of isogeometric analysis was the 
advent of Analysis-suitable T-splines (ASTS). ASTS are a mildly restricted sub- 
set of T-splines which are optimized to simultaneously meet the needs of design 
and analysis [EH E] ■ Linear independence of analysis-suitable T-spline blend- 
ing functions was established in 19J. An efficient local refinement algorithm 
for ASTS was developed in [3]. Later, it was shown that a dual basis, con- 
structed as in the tensor product settting, could be generalized to ASTS [2D] . 
This characteristic of ASTS is called dual compatibility. These results were then 
generalized to ASTS surfaces of arbitrary degree in |21) . 

In this paper we continue to develop the theory of ASTS spaces. Specifically, 
we provide a rigorous characterization of ASTS and show that the space of 
smooth parametric bicubic polynomials, defined over the extended T-mesh of 
an ASTS, is contained in the corresponding ASTS space. To accomplish this, 
the theory of perturbed ASTS spaces is developed and an ASTS dimension 
formula is established in terms of the topology of the extended T-mesh. We 
note that, unlike existing approaches, our dimension formula does not require 
that the T-mesh have any particular nesting structure. We then show that 
this characterization, coupled with the dual compatibility of ASTS, can be used 
to prove that ASTS spaces possess the same optimal approximation properties 
as tensor product B-spline spaces [22) . Next, we prove under what conditions 
two ASTS spaces are nested. This provides the theoretical justification for the 
analysis-suitable local refinement algorithm in [3 and provides a foundation 
upon which adaptive isogeometric analysis procedures may be developed in the 
future. 

This paper is organized as follows. Section [2] describes the T-mesh in index 
space, T-junctions, and the extended T-mesh. T-splines in the parametric do- 
main, blending functions, and T-spline spaces are defined in Section[3] Section[I] 
describes the conditions under which a T-spline is analysis-suitable. The theory 
of smoothly perturbed ASTS is developed in Section [5] Section [6] proves the 
conditions under which two ASTS spaces are nested. Using the characterization 
of ASTS spaces and dual compatibility several basic approximation results are 
proven in Section [7J Finally, Section [8] proves the dimension of ASTS spaces. 

2 The T-mesh 

An important object underlying T-spline spaces is the T-mesh. A T-mesh is used 
to determine T-spline basis functions and how they are arranged with respect 
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to one another. In other words, the mesh topology of the T-mesh determines 
the functional properties of the resulting space. In an attempt to adhere to a 
single notation and to reduce confusion, we define a T-mesh following much of 
the notation given in [THl 120] • For quick reference, [A] lists the most important 
notational conventions used throughout the text and where they are defined. 



2.1 Definition 

A T-mesh T is a rectangular partition of the index domain [to, to] x [n, n] , 
to, m, n, n € Z, where all rectangle corners (or vertices) have integer coordinates 
and all rectangles are open sets. Each vertex in T is a singleton subset of Z 2 . 
We denote all vertices of T by V. An edge of T is a segment between vertices of 
T that does not intersect any rectangle of T. We note that edges do not contain 
vertices and they are open at their endpoints. We denote all edges of T by E. 
Figure [l] shows an example of a T-mesh. The notation T 1 C T 2 will indicate 
that T 2 can be created by adding vertices and edges to T 1 . 
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Figure 1: A T-mesh. 



The valence of a vertex V <E V is the number of edges such that V is an 
cndpoint. We only allow valence three (called T-junctions) or four vertices. 
Note that valence two vertices, other than the four corners, are eliminated from 
the definition. 

The horizontal (resp., vertical) skeleton of a T-mesh is denoted by hS (resp., 
i>S), and is the union of all horizontal (resp., vertical) edges and all vertices. 
Finally, we denote the skeleton to be the union S = hS U vS. For a given vertex 
a = we define hS(a) := {k e Z : {k} x a C vS} and id (a) := {k G Z : 

a x {k} C hS}. We assume that these two sets are ordered. 



3 



We split the index domain R = [to, to] x [n, n] into an active region AR and a 
frame region FR such that R = FR U AR and AR = [to + 2, TO - 2] x [n + 2, n - 2] , 
and FR = Qto,to+ 2] U [to — 2, to]) x [n,n] U [to, to] x ([ri,n + 2] U [n — 2,n]). 
Note that both FR and AR are closed. 

A symbolic T-mesh |19) is created from a T-mesh T by assigning a symbol 
in Table [T] to each vertex in a tensor product mesh formed from the index 
coordinates, {to, . . . , to} x {n, . . . ,n} C Z 2 . The symbol is chosen to match 
the mesh topology of T. The symbolic T-mesh corresponding to the T-mesh in 
Figure [I] is shown in Figure [2] 



Table 1: Definition of possible symbols in a symbolic T-mesh 
Symbol Correspondence with T 

+ Valence 4 vertex, corner vertex, or valence 3 boundary vertex in T 

-, H, -L, T Oriented valence three vertex in T 

Vertical edge in T 
— Horizontal edge in T 

No corresponding vertex or edge in T 
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Figure 2: The symbolic representation of the T-mesh in Figure [T] 



2.2 Admissible T-meshes 

We say that a T-mesh is admissible if it satisfies three basic conditions. First, 
we require that S n FR contains the vertical segments {i} x [n, n] for i — to, to + 
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l,m + 2, to — 2, to — l,m and the horizontal segments [to, to] x {j} for j = 
Tki Tk + 1, 21 + 2, n — 2, n — 1, n. These horizontal and vertical lines are for basis 
function definition near the boundary. Second, we require that S fl AR contains 
the vertical segments {i} x [n,n] for i = to + 2, m + 3,to — 3, to — 2 and the 
horizontal segments [to, to] x {j} for j = n + 2, n + 3, n — 3, n — 2. Third, we 
require that for any two vertices V\ — {(«i,Ji)}, V% — {(*2, J2)} in V, such that 
Vi,V 2 C <9Q for some Q G T, if ii = i 2 (resp., ji = j 2 ), then {ii}x]ji, j 2 [c 
S (resp., ]«i,«2[x{ji} C S). From a practical point of view these are minor 
restrictions. The T-mesh in Figure [I] is admissible. 

We note that for convenience and simplicity, we often refer to only the active 
region of an admissible T-mesh when speaking of a T-mesh. In all cases, we 
assume that the frame region has an admissible topology. 

2.3 Anchors and T-junctions 

We define the anchors A(T) = {A e Vfl AR}. We denote the total number 
of anchors in T by n A . We define J C A(T) to be the set of all valence three 
vertices. These are called T-junctions. The symbols h, H, _L, T indicate the 
four possible orientations of a T-junction in a symbolic T-mesh. A T-junction 
Th <E J (resp., T v G J) of type h and H (resp., _L, T) and their extensions are 
called horizontal (resp., vertical). The solid white and red circles in Figure [l] 
are anchors and the red circles are T-junctions. 

2.4 Segments 

We define a segment to be a closed line segment of contiguous vertices and edges 
whose beginning and ending vertices are T-junctions (interior or boundary). 
Given two horizontal (resp., vertical) segments G\ defined over the intervals 
[«i, ji] x a and [12, J2] X b we say that < G\ if i\ < i 2 - We denote by hG 
(resp., vG) the collection of all horizontal (resp., vertical) segments, and by 
G = hG U vG the collection of all segments. We define hG(a) = hG n ([to, to] x a) 
and vG(a) — vG n (a x [n, n]). We assume these two sets are ordered. We 
denote the total number of segments in T by n G . We denote the total number 
of horizontal (resp., vertical) segments in T by n G (resp., n^ ). We denote the 
number of line segments in hG(a) (resp., vG(a)) by n^{a) (resp., n^(aj). 

2.5 The extended T-mesh 

T-junction extensions can be associated with each T-junction. For example, 
given a T-junction T = {(1, j)} e J of type h we extract from ftj({j}) four 
consecutive indices ii, . . . ,14 such that 1 = 13. We call ext e (T) — x {j} 

the face extension, ext? (T) =]z, ^4] x {j} the edge extension for such kind of 
T-junction. Similarly, we can define the face and edge extensions for the other 
kinds of T-junctions H, _L, T which are illustrated in Figure [3] 

We denote ext(T) = ext e (T) U ext^(T) the extension of T-junction T and 
the union of all horizontal (resp., vertical) face extensions by hext'(T) (resp., 
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Figure 3: The face extension (red) and edge extension (yellow) for four different 
kinds of T-junctions. 



vexv (T)), the union of all face extensions by earf/(T), and the union of all 
extensions (face and edge) by eart(T). We define the extended T-mesh, T ext , 
as the T-mesh created by adding to T all the T-junction extensions. In other 
words, Text — T U ext(T). We denote the total number of vertices in T ext by 
n ext rji^ gx^gjj^gj T-mesh corresponding to the T-mesh in Figure [l] is shown 
in Figure [4] 

Adding T-junction extensions to T may introduce three additional collections 
of vertices. The first, called crossing vertices and denoted by CV, is created from 
the intersection of crossing face extensions. In other words, 

CV = hext(J) Hvext{T). 

We denote the number of crossing vertices in T ext by n + . In Figure |H the 
crossing vertices are denoted by red stars. 

The second, called overlap vertices and denoted by OV, is created from the 
intersection of overlapping face extensions with S. In other words, 

OV = (( p| ext f (T h )) n vS) U (( p| ext f (T v )) n hS). 
T h ei t„sj 

We denote the number of overlap vertices in T ext by n~ . In Figure [4] the overlap 
vertices are denoted by green triangles. 
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Figure 4: The extended T-mesh corresponding to the T-mesh in Figure [T] The 
crossing vertices are denoted by red stars, the overlap vertices are denoted by 
green triangles, and the extended vertices are denoted by black squares. Notice 
that the inactive vertices (grey circles) are also regarded as extended vertices. 
The active vertices are denoted by hollow circles and the T-junctions are denoted 
by red circles. 



The third, called extended vertices and denoted by EV, is created from 
the intersection of face extensions and S while removing those vertices which 
already correspond to overlap vertices. Additionally, all non-anchor vertices are 
classified as extended vertices. In other words, 

EV = ((ext f (J) n S) \ OV) u (V \ A(T)). 

We denote the number of extended vertices in T ext by n* . In Figure [4] the 
extended vertices are denoted by black squares. 

3 The parametric domain and T-spline spaces 

Let H = (£ m , . . . , and II = (?y„, . . . , rfc) be two global knot vectors defined on 
the interval K. Interior knots may have a multiplicity of three while end knots 
may have a multiplicity of four. The global knot vectors define a full paramet- 
ric domain, fl C 1R 2 , where fl = [£ m ,fm] ® [r]n,Vn) and a reduced parametric 
domain, il C SI, where ft = [£m+3, £m-3]<8>[f7n+3, Vn-3]- The T-mesh in the para- 
metric domain is defined as the collection of non-empty elements of the form 
Q : Ci 2 [ x : %2 [ where Q =]ii,i2[x]ji,j2[G T. We denote those elements 
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where QC\ Cl ^ by Q. The extended T-mcsh in the parametric domain as well 
as all clement related concepts are defined similarly. Throughout this paper we 
use the index and parametric representation of a T-mesh interchangeably with 
the context making the use clear. 

For each anchor A = a x b E A(T) we define its horizontal (vertical) index 
vector hv{A) (vv(A), respectively) as a subset of hJ(b) (wJ(a), respectively) 
where hv{A) = . . . ,15) 6Z 5 contains five unique consecutive indices in hJ(b) 
with {13} = a. The vertical index vector, denoted by vv(A), is constructed in 
an analogous manner. We then associate a T-spline blending function N A (^r]) 
with anchor A. The T-splinc blending functions are given by 



n a (c,v) ■= B[E A ](OB\n A ](v) v(c, 77) e n 



(1) 



where B[E, A }(^) and B[U A ](r]) are the cubic B-spline basis functions associated 
with the local knot vectors 



"A — [Cii ! • • • > 615] G n 
U A = [r) h , . . . , r] h ] c n 



(2) 
(3) 



and hv(A) = ...,i 5 ) and vv(A) = (ji, . . . , j 5 ). 

Figure [5] illustrates the construction of a T-spline blending function cor- 
responding to anchor A = {(3,3)}. In this case, the local knot vectors are 

= [6,6, £3, £4,6,] and IF4 = [i7i, 772,»73,f?4,%]< 
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Figure 5: Inferring a T-spline blending function from a T-mcsh. 



A T-spline space T = T(T, S, II) is simply the span of the blending functions, 
N A , A e A(T). 



4 Analysis-suitable T-splines 



Analysis-suitable T-splines form a practically useful subset of T-splines. ASTS 
maintain the important mathematical properties of the NURBS basis while pro- 
viding an efficient and highly localized refinement capability. Several important 
properties of ASTS have been proven: 

• The blending functions are linearly independent for any choice of knots |19j . 



• The basis constitutes a partition of unity (see Corollary 6.7). 

• Each basis function is non-negative. 

• They can be generalized to arbitrary degree [2T] . 

• An afhne transformation of an analysis-suitable T-spline is obtained by 
applying the transformation to the control points. We refer to this as 
affine covariance. This implies that all "patch tests" (see are satisfied 
a priori. 

• They obey the convex hull property. 

• They can be locally refined [TJ [3] . 

• A dual basis can be constructed [501 HI] ■ 

Definition 4.1. An analysis- suitable T-spline is a T-spline whose T-mesh is 
analysis- suitable 119\/ . A T-mesh is said to be analysis-suitable if it is admissible 
and no horizontal T-junction extension intersects a vertical T-junction exten- 
sion. 

An analysis-suitable T-mesh is shown in Figure [6^,. The corresponding ex- 
tended T-mesh is shown in Figure |6]d. Notice that no horizontal extension 
intersects a vertical extension. The dual basis for an ASTS equips these spaces 
with a rich mathematical structure which we leverage in this paper (30]. 

Lemma 4.2. For a bicubic ASTS, each dual basis function, corresponding to a 
T-spline basis function Na{£,, rf) with local knot vectors = , . . . , £j 5 ] and 
n A = [r) h ,. . .,r] h ], is 

A.4 = A[£ 41 , . . . ,£j 5 ] ® •••) Vj 5 ] 

where A[£j x , . . . , & 5 ] and \[rjj 1 , ■ ■ ■ , r)j 5 ] are dual basis functions corresponding to 
univariate cubic B-splines WM whose knot vectors are , . . . , £j 5 ] and [rjj 1 , . . . ,rjj s 
respectively. 
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(a) Analysis-suitable T-mesh 



(b) The extended T-mesh 



Figure 6: An analysis-suitable T-mesh (a) and the extended T-mesh (b). 

5 Perturbed T-splines 

From a theoretical point of view, developing a complete and rigorous character- 
ization of T-spline spaces is complicated by the presence of zero knot intervals 
(especially near T-junctions) and overlap vertices. However, allowing both is 
important when T-splines are used as a tool in design and analysis. 

To overcome this difficulty, we develop the theory of the perturbed T-mesh 
(and resulting perturbed T-spline space). A perturbed T-spline can be used to 
prove properties about the original T-splinc. In other words, we will generate a 
perturbed T-mesh, establish the result in the perturbed setting, and then show 
that the result holds as the perturbation converges to the original T-splinc. 

5.1 Perturbed T-meshes 

A perturbed T-mesh is created by first generating perturbed global knot vectors, 
S[<5], II[<$], where 6 = (Si, . . . , 8 n ) T is a vector of perturbation parameters. A 
perturbed global knot vector is written as 

s[^ = (e[^,...,e^]„...,^k) 

where i = t(i,g) takes the index of the i th knot in S and the g th segment 
in vG({i}) and returns a unique index in the perturbed global knot vector. 
The knot values are initialized as £[8]iU,g) — £i- In other words, a knot index 
which corresponds to a vG({i}) which contains multiple segments in the T- 
mesh is repeated n^({i}) times. Notice that this operation induces an index 
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map fnr(i(i,g)) = i (resp., vtt) from the indices in the perturbed global knot 
vector onto the original global knot vector. The knot values are then perturbed 
using a small parameter S £ K as 

i=*+l 

where A£[6]j = c a{j) 6 = 6 a , if £ h7r(j ) - ^(j-i) = °, and is ec l ual to 6^0) ~ 
£hn(j-i)i otherwise. The constant, c a £ [0, oo). This same procedure is applied 
to n to form H[5]. The T-mesh, T, is then modified to form the perturbed T- 
mesh, T[S], by associating the vertices and edges contained in the g th segment of 
vG({i}) with knot £[<J]i(i,g)- Notice that the number of anchors does not change 
when forming T[S]. A perturbed T-spline space T[S] is a T-spline space formed 
from perturbed global knot vectors and T-mesh. A strictly perturbed T-mesh 
or T-spline space is one where c a £ (0, oo). 

A perturbation of an analysis-suitable T-mesh is shown in Figure [7| The 
analysis-suitable T-mesh is shown in Figure [7^, and the perturbed T-mesh is 
shown in Figure [TJa. Knot intervals are shown instead of knots for simplicity. 
Recall that a knot interval is simply the difference between adjacent knots in a 
global knot vector. Notice that the horizontal and vertical zero knot intervals 
have been replaced by non-zero knot intervals ai and 03. The vertical segments 
with T-junctions are perturbed resulting in a new knot interval a\. 
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(b) The perturbed T-mesh 



Figure 7: A perturbation of an analysis-suitable T-mesh T shown in (a) results 
in the perturbed T-mesh T[S] shown in (b). 



Proposition 5.1. IfT is analysis-suitable then T[S] is analysis-suitable. 

Proof. Suppose the extensions of two T-junctions ext(V\) = X {j 1 } (b 

or H) and extiy-z) = {i 2 } x (-L or T) in T[S] intersect. This implies 
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that jf < j 1 < jl and i\ < i 2 < i\. According to the construction of T[S], 
we don't change the order of the indices, so hn(jf) < hn^ 1 ) < /i7r(j|) and 
vn(i\) < vir(i 2 ) < vir(i\), i.e., there are intersecting extensions in T. □ 

Lemma 5.1. LetT be an analysis-suitable T-mesh. For every anchor, 

and horizontal index vector, hv({(i, j)})[8], in the perturbed T-mesh, T[S], hv({(hn(i),VTr(j))}) 
hir(hv({(i, j)})[5]) where hiv(hv({(i,j)})[6]) = ( hw(ii ),..., h/K{i§)). This also 
holds for the vertical index vectors. 

Proof. If the topological symbols corresponding to the indices {(ii,j), ■ ■ ■ , (is, j)} 
where ig € hv({(i, j)})[5], £ = 1, ... ,5, are not h or H then the result imme- 
diately follows since the indices are contained in a single horizontal segment. 
Thus, we only need to prove that the result holds when the topological symbol 
corresponding to {(«,j)} is h or H. Without loss of generality we assume it to 
be K Since, according to Lemma |5.1| the T-mesh, T[S], is analysis-suitable the 
symbols for the first two indices in hv({(i, j)})[5] can only be | or H. Thus, 
hv({(hn(i),vn(j))}) = hir(hv({(i, 3 )})[8}). ' □ 

Theorem 5.2. IfT is analysis- suitable andl~[S] is an offset perturbation then 
N A [S] -> N <A ) as5->0, Ae A(T[S}) where it {A) = {(/itt(«), im(j))}. 

Proof. Notice that if <5 — >• 0, £[6] z — > £,hir(i) for any i and i][S]j — > r] V7V (j) for any j. 



According to Lemma 5.1 the local knot vectors for each basis function iV^^] 
converge to those for N^i A \, A G A(T[<5]). Thus, according to ([Mi, Theorem 
4.36) we have the result. □ 



6 Refineability and nestedness 

We now explore the refineability and nesting behavior of analysis-suitable T- 
spline spaces. In other words, given two analysis-suitable T-splines spaces, T 1 
and T 2 , we establish the conditions under which T 1 C T 2 . We first establish 
basic refineability properties when the analysis-suitable T-mesh does not have 
any knot multiplicities or overlap vertices. Using the theory of perturbed T- 
splines, we then extend those results to encompass T-mcshcs which do have zero 
knot intervals and overlap vertices. 

Definition 6.1. The notation T 1 [S, T 2 ] denotes a perturbed T-mesh where T 1 C 
T 2 andT 1 [^,T 2 ] is created by removing those edges and vertices from the strictly 
perturbed T-mesh T 2 [8] which correspond to unperturbed edges and vertices in 
T 2 \ T 1 . By inspection, it is clear that T [6, T ], constructed in this way, is a 
perturbed T-mesh which satisfies Proposition \5. 1\ Lemma \5.1\ and Theorem \5 
and that T^tfjT 2 ] C T 2 [(5]. 

The construction of T 1 [^,T 2 ] is depicted in Figure [8] Two analysis-suitable 
T-meshes are shown in Figure [8^, and Figure [8p. Notice that T 1 C T 2 . The 
perturbed T-mesh T 1 [S, T 2 ] (shown in FigureT8p) is formed by removing the 
dotted lines (shown in Figure pp) from T 2 [ S] (shown in Figure pp) . 
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Figure 8: The construction of T 1 [^,T 2 ]. Two analysis-suitable T-meshes are 
shown in (a) and (b) such that T 1 C T 2 . The perturbed T-mesh T x [<$, T 2 ], 
shown in (c), is formed by removing the dotted lines from T 2 [(5], shown in (d). 



Definition 6.2. Given a T-mesh, T, with no knot multiplicities, and the cor- 
responding extended T-mesh, T ext , the homogeneous extended spline space is 
defined as 

s ext = {/ e c* 2 ' 2 (k 2 ) | f\ Qext £ ¥ 33 ,VQ ext c n, and/| R2 ^ = o} (4) 

where C 2,2 (R 2 ) is the space of bivariate functions which are C 2 -continuous in 
£ and r\ over all o/R 2 . P33 is the space of bicubic polynomials. The extended 
spline space is defined to be T ex t = Sext\n- 

Proposition 6.1. dimT ex t = dimS ex t 

Proof. We first prove that the dimension of S ex t is not less than the dimension 
of Text- Notice that for any function / £ S ex t, /|q £ Text- We now show that 
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the dimension of T ex t is not less than the dimension of S ex t- This is equivalent 
to showing that there is only one function in S ex t which is zero over f2. It is easy 
to see that the only function which is zero over ft must be zero over all of M 2 
since the minimum support of a cubic C 2 spline function is four intervals. □ 

Lemma 6.3. If the extended T-mesh, T ext , of an analysis-suitable T-mesh, T, 
has no knot multiplicities or overlap vertices, then T = T ex t- In other words, 
the analysis-suitable T-spline space, T, and the extended spline space, T ex t, are 
the same space. 

Proof. We have that T C T ex t (see [2D], Lemma 4.3), so the dimension of T is 



less than that of T ex t , which, according to Theorem |8.5[ Proposition 6T and 
Theorem |8.9| is the number of active vertices. Since the blending functions for 
anaysis-suitable T-splines are linearly independent the dimension of T is also 
the number of active vertices. Thus, the two spline spaces are identical. □ 

Lemma 6.4. Given two analysis-suitable T-meshes, T 1 and T 2 , neither of 
which has knot multiplicities or overlap vertices, ifT\ xt C T 2 xt , then T 1 C T 2 ■ 

Proof. Obviously, T e \ t C T^ xt . Since T 1 and T 2 are analysis-suitable, according 



to Lemma 6.3 T 1 C T 2 . □ 



Lemma 6.5. Given two analysis-suitable T-spline spaces, T 1 andT 2 , ifT 1 ^} C 
T*[8\, thenT 1 C T 2 . 

Proof. Suppose the perturbed T-spline space, T 1 ^], is spanned by the basis 
functions, JV^[d], and the perturbed T-spline space, T 2 [8], is spanned by the 
basis functions, iV|[($]. We have that 

N 1 a [S]=J2MNa[S})N 2 [8}. 

where Xb, B E A(T 2 [^]) are the dual functionals for the analysis-suitable T- 



spline basis Nb[S] as described in [5D]. According to Theorem 5.2 
hm N 2 [S] = N 2 (B) , lim N\ [8] = N\ A) . 

According to Theorem 4.41 in [23], \b{Na[8\) is bounded, so 

\im\ B {N A {8\)=4, 

8— >0 



i.e., ^ 



K (A) = J24n 2 



□ 



Theorem 6.6. Given two analysis-suitable T-meshes, T 1 andT 2 , if~[\ xt [8, T 2 ] C 
T 2 ext [S], thenT 1 CT 2 . 
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Proof. Obviously, 7^^. t [<5,T 2 ] C T 2 xt [5]. Since T^J^] has no knot m ultip licities 



or overlap vertices we have that T 2 [o~] = 7f xt [6] according to Lemma 6.3 Since 
the extended T-mesh, Tl xt [5, T 2 ], may have overlap vertices but no knot mul- 
tiplicities w e ha ve that T 1 ^, T 2 ] C 7^ t [5,T 2 ] according to Theorem 8.5 and 
Proposition 6.1 This immediately implies that T 1 C T 2 . □ 



Proposition 6.2. Every analysis-suitable T-spline space contains the space of 
bicubic polynomials. 



Proof. We have that any bicubic polynomial f[S] £ T[$\ according to Lemma 6.3 
Using the dual basis for T[S] we have that 

/[*! = EM/)^[«]. 

Thus, according to Theorem |5.2| and (|2"4"]. Theorem 4.41), as 6 — > 0, 

f = Y,c A N A . 

□ 

Corollary 6.7. Every analysis- suitable T-spline space forms a partition of 
unity. In other words, J^a ^a(^i]) = L V(£, rj) G fi. 



Proof. This immediately follows from Proposition 6.2 and the fact that A^(l) 



1. □ 



7 Approximation 

As described in [201 HI] approximation properties of analysis-suitable T-splines 



are directly linked to Proposition 6.2 In other words, having the bicubic poly- 
nomials in the T-spline space is the minimal requirement to obtain an 0(h 4 ) 
convergence rate in the mesh size. 

Following the approach in 22, 23 121], the dual basis for an analysis-suitable 
T-spline space, T, can be used to construct a projection operator, P : L 2 (fi) — > 
T, where 

nm,v)= E *AV)N A (£,r,) VfeL 2 (Cl),V(Z,v)e&- 

AeA(T) 

We denote the open support of a T-spline basis function by Qa C fi, and the 
extended support of an element Q by fig C fi, where 

fi Q = IJ Qa, A(Q) = {A G A(T) : Q A n Q ± 0}. 

AGA(Q) 

We will denote by -R(fig) the smallest rectangle in fi containing fig and . 
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Proposition 7.1. Given an analysis-suitable T-spline space, T, the projection 
operator P is (locally) h-uniformly continuous in the I? norm. In other words, 
there exists a constant C independent of T, S, II such that 

WHDWlhq) < C\\f\\mn Q ) VQ g T ext ,Mf G L 2 (Q) 

Note that the constant C may depend on the polynomial degree. 

Proof. The result follows immediately from ([2Q], Proposition 5.4) and Propo- 
sition lOl □ 



Proposition 7.2. Given an analysis- suitable T-spline space, T , there exists a 
constant C independent o/T, S,II such that for r £ [0,4] 

11/ -nDW^Q) < C"(h mo) ) r \f\ Hr{R{no)) VQ G T ext ,Vf € L 2 (Q) 

where h^Q^ denotes the diameter of R(Qq). Note that the constant C may 
depend on the polynomial degree. 

Proof. The result follows immediately from ([50], Proposition 5.4) and Propo- 
sition lO □ 



8 Dimension 



In this section, we develop a dimension formula for polynomial spline spaces 
defined over the extended T-mesh in the parametric domain of a T-spline and 
establish the connection between this dimension formula and analysis-suitable 
T-spline spaces. The dimension formula, written only in terms of topological 
quantities of the original T-mesh, is an essential ingredient in establishing the 
refineability properties in Section [6] and the approximation results in Sectio n [7 
for analysis-suitable T-splines. The essential results are proven in Theorems |8. 5 
and EH 

Unlike existing approaches, our dimension formula does not require that the 
T-mesh have any nesting structure. Of critical importance is how this dimension 
formula can be directly related to analysis-suitable T-spline spaces which can 
then be used to construct a simple set of basis functions for the spline space 
which are compatible with commercial CAD and analysis frameworks. 



8.1 Smoothing cofactor-conformality method 

We use the smoothing cofactor-conformality method [551 HI] to transform the 
smoothness properties of S ext into a linear constraint matrix, M. This con- 
straint matrix is then analyzed to determine the dimension of S ex t- We recall 
that the spline space, S ex t, is defined using the extended T-mesh, T ext , corre- 
sponding to a T-mesh, T, which does not have any knot multiplicities. 
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8.1.1 Vertex and edge cofactors 



As shown in Figure [9j for any vertex, Vij — (£i,T]j) G Text, the surrounding 
bicubic polynomial patches are labeled, Pij(£,r]), k — 0,1,2,3. If the vertex, 
Vij, is a T-junction, then Pij(£, rj) = p^j 1 ^, rj) for some k. Since Pij{£, v) an d 
Pij{€,rj) are C 2 -continuous there exists a cubic polynomial 7^(17), called the 
edge cofactor, such that 

pUt,v)-P° i , j (t;,v)=iUm-{;i) s - (5) 



7- • 



7i, 



Mi. 7 



On 



7 2 ' 



7 2 - 



— ^— Q — 



7 2 - 



-Or 



(a) 



(b) 



(c) 



Figure 9: The smoothing cofactors around a vertex. 

Similarly, there exists cubic polynomials, jljii]), /J-jjiOi an d KjiOi sucn 
that 

pUZ, v) - pUZ, v) = rijWn ~ ^) 3 ' 

pL-(£, »?) - »?) - -T^toX* - &) 3 , (7) 

Pijfav) -PU^V) = -rf,j(0(ri-Vj) 3 - (8) 

We note that if two patches are identical the edge cofactor is zero. Combining ([5| 
- Q gives 

hUv) - lUvMt - &) 3 = - ^-(O)^ - mf- (9) 

Since (£ — £i) 3 and (77 — 77^ ) 3 are prime to each other there exists a constant, 
dij, called the vertex cofactor, such that 

lh fa) ~ lh fa) = dij fa - V,) 3 , lA d (0 - Mi.j (0 = dij « ~ 6) 3 - (10) 
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Figure 10: The smoothing cofactors along a horizonal edge segment. 



8.1.2 Assembling the constraint matrix, M 



Referring to Figure 10 consider a horizontal segment G^ with k + 1 vertices 



and k edge cofactorsTTJsing ( 10 ) we have that 



(11) 

(12) 



o-/4,,- = ^,j(£-&J 3 , 



and 



U+1,3 



H,3 ' 



i = l,...,k. 



Summing (11) - (13 1 and using (14) results in the linear system 



(13) 
(14) 

(15) 



We call the solution space, denoted by M^[G^], for this linear system the edge 
conformality space. Similarly, for a vertical segment G\ we have that 



(16) 



where the solution space is denoted by W[G|]. By (15) and (16), one immedi- 
ately has that 

Lemma 8.1. // each £j £ and r]j e are different, then the dimension of W[Gj] 
and W[G!-] are k — 3 and I — 3 respectively. 



The linear systems ( 15 ) and ( 16 1, associated with the horizontal and vertical 



segments in T ex t, can be assembled into the global system 

MD = 



(17) 



where D = [di, d2, ■ ■ ■ , d n ^t ] T is a column vector of all vertex cofactors in 
T ext and M is a 4n seg x n ext real matrix. Each edge conformality condition 
corresponds to a submatrix consisting of 4 rows of M and each vertex cofactor 
corresponds to a column of M. 
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Lemma 8.2. The dimension of S ex t is the nullity o/M, i.e., the dimension is 
n ext minus the rank o/M. 

Proof. Since the continuity constraints in S ex t have been converted into the lin- 



ear system in (17) using the smoothing cofactor-conformality method (see |25j). 
the dimension of S ex t is the dimension of the null space of M, i.e., the dimension 
is n ext minus the rank of M. □ 



8.2 Simplifying the constraint matrix, M, and T ext 

It is possible to simplify the constraint matrix, M, and the topology of the ex- 
tended T-mesh in the parametric domain, T ext , such that the null space of M 
is undisturbed. To remove a vertex from T ext means we delete the correspond- 
ing column from M and to remove a segment from T ext means we delete the 
appropriate submatrix from M. We form the reduced constraint matrix M by 
removing the eight segments and contained vertices [^miCm] x {Vn}, [£mj£m] x 

{Vti+i}, Km, Cm] X {rfc-l}, [f m j 6n] X {Vn}, {Cm} X [n^Vn], Um+l} X [??„ , 7fa] , 

{£m-i} x [rjn, rjn], and {£,m} x [rjn, rjn}. We denote the T-mesh after the removals 
by Text and th e nu mber of vertices and segments in T ext by n ext and n G , respec- 
tively. Figure 11 shows the simplified extended T-mesh T ext for the extended 
T-mesh in Figure [4] The vertices and segments which remain after the removal 
process have corresponding entries in M. 
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Figure 11: The simplified extended T-mesh T ext for the extended T-mesh in 
Figure |4j 



Lemma 8.3. The dimension of the null space o/M is the same as that for M. 
Proof. The vertex cofactors which correspond to the removed corner vertices 



can be uniquely determined by applying (15) to the four horizontal removed 
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segments or by applying ( 16 ) to the four vertical removed segments. To establish 



the result we need to show that the constraints corresponding to the four vertical 
removed segments can be derived from the constraints corresponding to the four 
horizontal removed segments. We have that 



°=E 

m— C 
I 

= E 



k 



m=0 



(£-6 J 3 



E 

1=0,1,1-1,1 



.m=0 



(18) 
(19) 
(20) 



Equation ( 18 1 involves the sum of all edge conformality conditions for the 



horizonal edge segments. 



Equation (20 1 holds because the linear systems for 



the other vertical segments are satisfied. Since (77 — r]j e ) , I — 0, 1,1 — 1,1, 
form a basis for a linear space of polynomials with degree less than four, 
12m=o di m ,ji(£ ~ £i m ) 3 = 0, for m = 0, . . . , fc. In other words, the constraints 
for the four vertical removed segments can be derived from the other con- 
straints. □ 

Remark 8.4. There are many possible simplification techniques which could 
have been used. This simplication technique was chosen because it leaves the 
null space of M undisturbed and each resulting segment in M contains exaclty 
four extended vertices. 



Theorem 8.5. I/M has full column rank, then the dimension of S ex t is 

dim^e^t = n a + n + + n~ 



(21) 



where n a is the number of active vertices in T and n + and n are the number 
of crossing and overlap vertices, respectively, in T ext . 

Proof. Since there are n ext and n G vertices and segments, respectively, in T ext , 
M is a 4n G x n ext matrix. Since M has full column rank the dimension of S ex t 



is n ext 



4n . As every segment in T ext has exactly four extended vertices and 



these four extended vertices are not extended vertices for any other segment, 
the number of extended vertices in T ext is An G . Thus, 



dimSe 



n + n + + n 



(22) 



□ 



8.3 Rank of the constraint matrix M 

Since every extended vertex in T ext is an extended vertex in exactly one segment, 
the matrix M has more columns than rows, i.e., n ext > 4n G . After arranging 
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the order of edge conformality conditions and the order of vertex cofactors, an 
appropriate partition of the linear system of constraints, M D = 0, is 

= (23) 

where Mi is a 4n G x 4n G matrix and Mi is a 4rT G x (n ext — 4n G ) matrix, Di 
is a vector of the first 4n G vertex cofactors, and D2 is a vector of the remaining 
vertex cofactors. 

Definition 8.6. A simplified extended T-mesh is called diagonalizable if we can 
arrange all the segments Gj,i = 1, . . . , n G such that the number of vertices on 
segment Gj but not on segment Gi,i < j is at least 4- 

Lemma 8.7. If a simplified extended T-mesh is diagonalizable, then the matrix 
M has full column rank. 

Proof. Since the simplified extended T-mesh is diagonalizable, the segments 
can be ordered as described in Definition |8.6| Given this ordering of segments, 
for i — 1, . . . ,n G , we place the edge conformality conditions corresponding to 
segment Gi in rows 4(rT G — i) + 1 through 4(n G — i) + 4 of M, and place the 
first four vertex cofactors which appear in Gi but not in Gj, j < i in columns 
4(n G — i) + 1 through 4(n G — i) + 4 of Di. Then the matrix Mi is in upper 
block triangular form and according to Lemma |8.1| each diagonal block 4x4 
matrix is full rank, thus matrix Mi is obviously of full rank. □ 

Lemma 8.8. T ext is diagonalizable if and only if for any set of segments there 
exists at least one segment in the set that has at least four vertices which are 
not in the other segments in the set. 

Proof. Assume the T-mesh is diagonalizable but there exists a set of segments 
{Gij ,j = 1, ■ • ■ , s} such that any segment in the set has at most three vertices 
which arc not on the other segments in the set. Without loss of generality, 
assume the diagonalizable segment ordering for T ext is Gi,i — 1, . . . 1 n G . Let 
A; be the maximal index for all ij,j — 1, . . . ,s and consider the set {G;,i — 
1, . . . , fc}. Since G& has at most three vertices which are not in the segments 
Gj. , j — 1, . . . , s, there exists an index i,i < k such that the number of vertices 
on segment Gk but not on segment G; is at most three. This violates the 
assumption that T ext is diagonalizable. 

Suppose for any set of segments in T ext there exists at least one segment 
which has at least four vertices which are not on the other segments in the set. 
For the set containing every segment, according to the assumption, there exists 
one segment, G s g, which has at least four vertices which are not on the other 
segments. Now, removing G s g from the set, we have that in the set of remaining 
segments there exists one segment, G h g_ 1; which has at least four vertices which 
are not on the other segments. Continuing this process, we can arrange all the 
segments Gi,i — 1, 2, . . . ,n G such that the number of vertices on segment Gj 
but not on segment Gi, i < j is at least 4. Thus T ext is diagonalizable. □ 



[ Mi I M 2 ] 
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Theorem 8.9. For an analysis- suitable T-mesh, matrix M has full column 
rank. 

Proof. For an analysis-suitable T-mesh we will prove that the corresponding 



simplified extended T-mesh is diagonalizable. Otherwise, according to Lemma 8.8 
there exists a set of segments such that each segment in the set has at most three 
vertices which are not on the other segments in the set. It is evident that the 
set must contain horizonal segments. Otherwise, any vertical segment violates 
the assumption because it must have at least four vertices which are not in the 
other segments in the set. Let Gi be the bottommost horizonal segment in the 
set (if there is more than one such segment choose one of them). Since Gi has 
at most three vertices which are not on the other segments in the set, one of the 
four extended vertices of Gi must lie on a vertical segment, Gj. Now, referring 
to Figure |12[ since the T-mesh is analysis-suitable, Gj must have two anchor 
vertices whose vertical index coordinate is less than a. Otherwise, the T-mesh 
is not analysis-suitable due to intersecting T-juction extensions. Additionally, 
Gj must have two extended vertices whose vertical index coordinate is less than 
a. Thus, there are four vertices in Gj which do not belong to any other segment 
in the set which contradicts the assumption. □ 



-o Gi = [h,k 2 ] x {a} 



G j = {b}x[£ 1 ,£ 2 ] 



Figure 12: A schematic for Theorem |8. 9 



9 Conclusion 

We have established several important properties of analysis-suitable T-splines. 
We developed a characterization of ASTS spaces and proved that the space of 
bicubic polynomials, defined over the extended T-mesh of an ASTS, is contained 
in the corresponding ASTS space. We then proved the conditions under which 
two ASTS spaces are nested. This provides the theoretical foundation for the 
analysis-suitable local refinement algorithm in [3j. Using the characterization 
of ASTS we then proved several basic approximation results. Additionally, we 
developed the theory of perturbed ASTS and a simple mesh-based dimension 
formula which is written in terms of the vertices in the extended T-mesh of a T- 
spline. Both of these developments were critical for the proofs in this paper and 
may have important applications in other contexts. While the developments in 
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this paper are restricted to bicubic surfaces the extension to arbitrary degree 
should be straightforward. 
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A Notational conventions 

In Tables [2] and |3] the notational conventions used throughout the text are listed 
as well as the section where they are defined. 
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Table 2: Notational conventions used throughout the text and where they are 
defined. 



Symbols 


Description 


Section 


T 
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2.1 




V 


All vertices in a T-mesh 


2.1 




E 


All edges in a T-mesh 


2.1 




S 


The skeleton of a T-mesh 


2.1 




hS OS) 
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2.1 




Q 


An element of a T-mesh 


2.1 




G 


All segments in a T-mesh 


2.4 




n G 
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hG (uG) 
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R 
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AR 


The active region of an index domain 
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FR 
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~T ext 


An extended T-mesh 


2.5 
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earf(T) 
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ext(T) 
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hextf(T) vext^T 


All horizontal (vertical) faces extensions in a T-mesh 


2.5 




CV 


All crossing vertices in T ext 


2.5 




n+ 


The number of crossing vertices in T ext 


2.5 




OV 


All overlap vertices in T ext 


2.5 




n~ 


The number of overlap vertices in T ext 


2.5 




EV 


All extended vertices in T ext 


2.5 




n* 


The number of overlap vertices in T ext 


2.5 




n ext 


The number of vertices in T ext 


2.5 
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Section 



A(T) 

n A 
J 
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(n A ) 

n 

Q 
n 

Q 

R^q) 

Qa 

r 

Sext 
Text 



All anchors in a T-mesh 
The number of anchors in a T-mesh 
The set of T-junctions in a T-mesh 
The horizontal (vertical) local index vector for N A (^,r]) 
The horizontal (vertical) local knot vector for Na(£,,t]) 5 
The full parametric domain of T |3 
A T-spline element in the full parametric domain 
The reduced parametric domain of T 
A T-spline element in the reduced parametric domain 
The extended support of Q 
The smallest rectangle containg fthatQ 7 
The T-spline blending function for anchor A 3 
The support of N A 7 
The dual basis function associated with N A 3 
A T-spline space 3 
A homogeneous extended spline space |5 
An extended spline space 



s[($](II[<5]) A horizontal (vertical) perturbed global knot vector 

i = i(i,g) (j = i(j,g)) the index of the g th segment in hG({i}) (vG({j})) 



T(S) 
T{8) 
T(5,T* 



the map from perturbed index i{i,g) to index i 

A perturbed T-mesh 

A perturbed T-spline space 

A perturbed T-mesh generated from T l (5) 



5 
5 

I 

5 
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